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Abstract-Inordinate localized temperature rise in the power 
transformer causes subsequent thermal breakdown. To prescribe 
the limits of short-term and long-term loading capability of a 
transformer, it is necessary to estimate the hottest spot 
temperature (HST) of transformer. This paper proposes the 
steady state temperature distribution of the power transformer 
windings. Oil in the transformer is assumed nearly 
incompressible and oil properties such as thermal conductivity, 
special heat, viscosity, and density vary with temperature. Finite 
difference method is used for numerical solution. The selected 
model for simulation is a 50KVA, 20 kV/400V oil natural, and air 
natural cooling (ONAN) power transformer. A Comparison of 
the author's results with those obtained from finite integral 
transform and experimental test confirms the validity and 
accuracy of the proposed method. 
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I. INTRODUCTION 

Thermal impact is a limiting factor for the transformer 
operation [1], [2]. Therefore, the knowledge of the 
temperature, especially the hottest spot temperature, is of high 
interest. If the temperature rise goes beyond the permissible 
value, in order to preserve the insulation from deterioration, 
the load of transformer must be reduced or an auxiliary 
transformer is used. For an oil-immersed transformer, the oil 
surrounds the transformer body. Oil is a nearly incompressible 
fluid and density changes due to temperature rise, therefore oil 
moves in the transformer. The heat transferred by convection 
is the most important method of heat transfer. The analytical 
solution of convection equation is normally difficult and 
sometimes it is impossible due to the complexity of the 
geometry. 

The basic criterion for transformer loading is the hottest 
spot temperature of the solid insulation. It must not exceed the 
prescribed value in order to avoid insulation faults. A 
procedure of hottest spot temperature calculation is given in 
the international standards [3]-[5]. The algorithm for 
calculating the hottest spot temperature of a directly loaded 
transformer using data obtained in a short circuit heating test 
is given in [6], [7]. The usage of average heat transfer 
coefficient is typical in a transformer designing process to 
calculate needed number (area) of cooling surfaces. 

In this paper, a procedure for obtaining the temperature 
distribution in the power transformer is proposed. For this 
reason energy and Navier-Stokes equations are solved using 
finite difference method. Therefore, a code has been provided 
under MATLAB software. The model can be used for 
temperature calculation on the arbitrary change of current and 
outside air temperature. In the paper, thermal model is in 



section II. Results and discussion of the proposed work have 
been provided in Section III. 

II. THERMAL MODEL 

Fig. 1 shows cross section of a 50KVA, 20kV7400V power 
transformer in two dimensions. Dimensions and specifications 
of the power transformer have been summarized in TABLE I. 
TABLE II shows the losses. Energy equation for Newtonian 
incompressible fluid such as oil in two dimensions is [8] 
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Where 

T: temperature 

V v : velocity in x direction 

V y : velocity in y direction 

k: thermal conductivity 

p : density 

C p : special heat 

ju: viscosity 

In (1), oil properties vary with temperature [9], [10]. The 
temperature dependence of oil properties are given: 
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Where fi is volumetric expansion coefficient. The nine 
constants for transformer oil have been listed in TABLE III. It 
is generally valid for all transformer oils that the variation of 
the oil viscosity with temperature is much higher than the 
variation of other oil properties [10]. Thus, all oil physical 
properties except the viscosity can be replaced by a constant. 
However, in this paper has been considered the influence of all 
oil properties. In (1) velocity is unknown, and then we must 
solve Navier-Stokes equations. Navier-Stokes equations in 
two dimensions for incompressible fluid are [8] 
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Continuity Equation: 

S^ + ^ = (7) 

dx By 

Momentum Equations: 
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Where P is pressure and g is gravitational acceleration. 



Table III. Oil constants 



Figure. 1 Three phase power transformer in two dimensions 
Table I. Dimensions and specifications of the proposed Power 

TRANSFORMER 



Core diameter (cm) 


10 


Width of each window 


12 


(cm) 




Height of window (cm) 


28 


Thickness of LV winding 


1.2 


(cm) 




Thickness of HV winding 


2 


(cm) 




Rated power 


50KVA 


HV voltage 


20kV 


LV voltage 


400V 


HV current 


1.44A 


LV current 


72A 


Table II .Power transformer losses 


Losses (w) 


Value 


Core 


158 


DC of LV windings 


384 


DC of HV windings 


534 



Eddy currents of LV windings 
Eddy currents of HV windings 



64 



Oil constant 


Transformer oil 


a. 


0.0000013573 


a 2 


2797.3 


a 3 


1960 


(■U 


4.005 


as 


887 


a 6 


-0.659 


a? 


0.124 


a s 


-0.0001525 


a 9 


0.00086 



The heat conduction equation for core and windings is 
written 



d 2 T d 2j 

k x — r+K — r + 2 : 

dx By 



(10) 



Where k x and k y are thermal conductivity in x and y 
directions respectively. The term Q is the volumetric heat 
source function and has been modified here to take care of 
variation of electrical resistance of copper with temperature. 
The heat source term Q can be of the form [9], [10] 



Q=Q Q -H + a c -(T-T ref )] 



(11) 



Where a c is the temperature coefficient of electrical 
resistance of copper wire. With this representation, the 
function Q becomes temperature dependent, distributed heat 
source. Thermal conductivities are unequal in different 
directions. Thermal conductivity has been treated as a vector 
quantity, having components in both radial and axial direction. 
Resultant thermal conductivity of the system is [9] 



K = Jk! + k 2 , 



(12) 



Where 
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Term K represents resultant thermal conductivity of 
insulation and conductor system. In this paper, the bottom oil 
rise over ambient temperature has been calculated as 



,i;r, +i 
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(13) 



Where 6 U is the bottom oil temperature rise over ambient 
temperature, 9fi is the full load bottom oil temperature rise 
over ambient temperature obtained from an off-line test, Ri is 
the ratio of load loss at rated load to no-load loss. The variable 
is the ratio of the specified load to rated load. 
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/. = 



rated (14) 

The exponent n depends upon the cooling state. The 
loading guide recommends the use of n=0.8 for natural 
convection and n=0.9-1.0 for forced cooling. 

m. RESULTS AND DISCUSSION 

A software program has been provided for finite difference 
solution of energy and Navier-Stokes equations using 
MATLAB software. Depth cross section of the proposed 
power transformer has been shown in Fig.2. Mesh size used 
for simulation is about 150x300 (150 in horizontal direction 
and 300 in vertical direction). We assume inlet oil temperature 
is temperature reference. Therefore, we have 
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Where T hoi i is the bottom oil temperature. Fig. 3 shows 
temperature distribution from LV winding in one per unit (p.u.) 
load with oil natural cooling (ON). It can be pointed out that 
for LV winding, maximum temperature location is around 
80% of winding height from the bottom and at about 50% of 
radial thickness of the layer. Temperature distribution from 
HV winding has been shown in Fig.4. It can be observed that 
the maximum temperature occurs in the neighborhood of 55% 
of the axial and 50% of the radial thickness of the layer. Fig. 5 
shows HST rise over bottom oil temperature versus load at 
T am b=25( C). Table IV shows comparison of the proposed 
method with finite integral transform used in [9] and 
experimental measurement in [11]. 



Figure. 2 Depth cross section of the proposed power transformer 




Height of LV winding(m| 



Radial thickness ot LV winding(m) 




Radial thickness of HV windingjm) 



Figure. 4 Temperature distribution (1 p.u.) of HV winding at Tamb=25(°C) 




Figure. 5 HST Rise over bottom oil temperature versus load at Tamb=25(°c) 



TABLE IV. HST (°C) MAGNITUDES AT T A „ B =25(°C) 



Load (p.u.) 


Proposed 


Analytical [9] 


Experimental [11] 


0.6 


51.72 


53 


50 


0.7 


58.35 


56 




0.8 


65.71 


66 




0.9 


73.78 


75 




1.0 


82.54 


80 


81 


1.1 


91.99 


93 


90 


1.2 


102.15 


101 


99 


1.3 


113.01 


115 




1.4 


124.59 


125 




1.5 


136.91 


139 





IV. CONCLUSION 

In this paper, we solved the heat transfer partial differential 
equations to find temperature distribution in power 
transformer windings. The results seem to correspond 
reasonably well with results of calculations and actual tests [9], 
[11]. The authors wish to point out that the IEEE loading 
guide offers relations for the calculation of the HST based on 
per-unit load. The formulations tend to ignore the possibilities 
of two transformers that are rating identical but have a 
different winding structure and varying heat loss/unit volume. 
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Figure. 3 Temperature distribution (1 p.u.) of LV winding at Tamb=25(°C) 
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